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Abstract: We study the quark number susceptibility in holographic QCD with a finite 
chemical potential or under an external magnetic field at finite temperature. We first 
consider the quark number susceptibility with the chemical potential. We observe that 
approaching Tc from high temperature regime, Xg/^^ develops a peak as we increase the 
chemical potential, which confirms recent lattice QCD results. We discuss this behavior in 
connection with the existence of the critical end point in the QCD phase diagram. We also 
consider the quark number susceptibility under the external magnetic field. We predict 
that the quark number susceptibility exhibits a blow-up behavior at low temperature as 
we raise the value of the magnetic field. We finally spell out some limitations of our study. 
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1. Introduction 

Understanding the QCD phase diagram is one of the important problems in modern the- 
oretical physics; see ||l[ for some recent reviews. An important feature of the QCD phase 
diagram is the existence of the critical end point (CEP) where the first order phase tran- 
sition, from the hadronic phase to the quark gluon plasma (QGP), ends. Based on the 
fact that at the critical point thermodynamic observables such as susceptibilities diverge 
and the order parameter fluctuates on long wavelengths, basic ideas to observe CEP are 
suggested in Q. 

One of the important signals of QGP formation is the fluctuation of conserved charges 
such as baryon number or electric charge ^. The quark (or baryon) number susceptibil- 
ity, which measures the response of QCD to a change of the quark chemical potential is one 
of them |3|, |5|. There have been many studies to calculate the quark number susceptibility 
in various model studies ^ and lattice simulations |l^, |ll|, |12|. In [13|, the quark 



number susceptibility at finite temperature is studied with AdS/QCD models. The exis- 
tence of a peak in the quark number susceptibility near Tc is confirmed by recent lattice 
QCD calculations based on the Taylor expansion with respect to the quark (or baryon) 



chemical potential |14, 15]. This implies the existence of CEP, at which the first order 



phase transition terminates in the (//g,r) plane of the QCD phase diagram, see [16| for a 
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review. Lattice QCD estimates that the location of the CEP would be = 165—175 MeV 
and /i^ = 250 — 400 MeV |16]. The existence of the CEP has been also investigated in 



various QCD-based model studies |17] 



The AdS/CFT correspondence [18, |l^, ^ is a powerful tool to investigate strongly 
coupled gauge theories including critical phenomena. Using this correspondence, we can 
obtain physical quantities in gauge theories from calculations in gravity side. Although 
the correspondence between QCD and gravity theory is not known, we can obtain much 
insights by using this correspondence. In fact, many properties are shared by various gauge 
theories, especially in hydrodynamic limit. Hydrodynamic properties can be read off by 



introducing a small perturbations into the bulk fields pi|, 22]. This gives small fluctuations 



to the fluid in the gauge theory side. Physical quantities, like transport coefficients, can 



be calculated by considering the responses to these small perturbations; see |23, 24] and 
references therein. 

In order to study dynamics of quarks (or baryons) in the gauge theory side, we have to 
introduce the U{1) baryon symmetry in the gauge theory side. According to the Gubser- 
Klebanov-Polyakov/Witten relation 20|, a global symmetry in the gauge theory side 



corresponds to a gauge symmetry in the gravity side. To consider the ?7(1) baryon sym- 
metry, we have to introduce a U{\) gauge field in the bulk. This kind of models can be 
constructed from the string theory setup. In D3/D7 setup, for example, D7-branes are 
considered as the flavor brane and gauge fields on D7-branes play the role of the bulk 
gauge field. This model has N = 2 supersymmetry, and hence, the corresponding gauge 
theory \s N = 2 super symmetric QCD (SQCD). Even though the real QCD does not have 
supersymmetry, this discrepancy is expected to be solved by universality mentioned above. 

In this paper, we study one of the promising QGP probes by using the AdS/CFT 
correspondence. The primary goal is to calculate the quark number susceptibility at finite 
temperature with a finite quark chemical potential. Our approach has some ups and 
downs. Contrary to the lattice QCD considerations, we can calculate the susceptibility 
with arbitrary values of the chemical potential. However, we are not able to study the 
quark number susceptibility in confined phase, which will be discussed at the end of the 



section 3.1. Moreover, our study based on AdS/CFT is inherently suffering from \/Nc 
corrections. Unfortunately, a systematic way of collecting all those corrections has not been 
established. With this caution in mind, we first revisit the quark number susceptibility at 
finite temperature, and then we generalize it with the chemical potential. For this we work 
in the AdS black hole and Reissner-Nordstrom-AdS backgrounds. 

In addition to this, we evaluate the quark number susceptibility at finite temperature 
with a constant magnetic field. The basic motivation is due to the observation that the 
constant magnetic field enhances the dynamical chiral symmetry breaking ]|2^. On top 
of it, recently it is argued that sufficiently large magnetic fields are created in heavy ion 
collisions 1^6|, and so our study may be tested in a terrestrial experiment. Therefore it is 
interesting to study the phase diagram of QCD as a function of the magnetic field, and so 
the quark number susceptibility. 
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2. Retarded Green functions and quark number susceptibility 



In this section, we briefly summarize the prescription for the Minkowskian correlator in 
the AdS/CFT correspondence, and then define the quark number susceptibility through 



the correlator. We here follow the prescription proposed in |21]. We work on the following 
5D background, 

ds^ = g^udxf^dx" + guuiduf, (2.1) 

where and u are the 4D and radial coordinate, respectively. We refer the boundary at 
u = and the horizon at ti = 1. Let us consider a solution of an equation of motion in this 
5D background. Suppose the solution is given by 

<t>{u.x)=j^^e^'^h{u)<tP{k), (2.2) 

where fk{u) is normalized such that /jt(0) = 1 at the boundary. After putting the equation 
of motion back into the action, the on-shell action might be reduced to surface terms 



u=l 

(2.3) 

u=0 



Here, the function Q{k,u) can be written in terms of f±k{u) and duf±k{u)- Accommodat- 
ing Gubser-Klebanov-Polyakov/Witten relation ^ to Minkowski spacetime, Son and 
Starinets proposed the formula to get the retarded Green functions. 



G^{k) = 2g{k,u) 



(2.4) 

u=0 



where the incoming boundary condition at the horizon is imposed. In this paper, we 
consider correlators of U{1) currents J^{x), where J^{x) is the vector current of quark field 
or quark number current. Now we define the precise form of the retarded Green functions 
which we discuss later: 

,{k) = -ij d^x e-'''^9{t){[J^{x), MO)]). (2.5) 

The quark number susceptibility was proposed as a probe of the QCD chiral phase 
transition at zero chemical potential |^ |5|, 

--IS- <-) 

Later it has been shown that the quark number susceptibility can be rewritten in terms of 
the retarded Green function through the fluctuation-dissipation theorem |^, 

XqiT, i^) = - lim ReiCt t{uj = 0, k)) , (2.7) 

where y{oj,k) is the retarded Green function, which is defined in (p.5|). 
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3. Reissner-Nordstrom-AdS background 

In this section, we briefly review the Reissner-Nordstrom-AdS (RN-AdS) background in 5D 
spacetime first. Afterwards, we explain the interpretation of RN-AdS spacetime in terms 
of the D3/D7-brane setting. 

We consider the Einstein-Maxwell action with negative cosmological constant: 

S=^J d^x^(^R - 2a) - d^xV=5^mn^", (3.1) 

where we denote the gravitation constant and the cosmological constant as = SvrGs and 
A, respectively. The U{1) gauge field strength is given by Fmn{x) = dmAn{x) — dnAm{x). 
The RN-AdS background is a solution of the equations of motion derived from this action, 
and it is given by 

d^' = S f - /(0(di)' + + -^Mr)\ (3.2a) 



P \ J r^f{r 

At = -^ + 11, (3.2b) 

with 

/W = l- — + — ' A = -^, e=^K. 

The parameters /, m, and q are the AdS radius, mass and charge, while Q and fi are 
C/(l) charge and constant which is interpreted as the chemical potential. Since the gauge 
potential At{x) must vanish at the horizon, the charge Q and the chemical potential n are 
related as ^ 

H = -2- (3.3) 

' + 

The horizons of the RN-AdS black hole are located at the zero for /(r)^, 

fir) = ^-^ + '^ = ^,{r'-rl){r'- r^) (r^ - r^) , (3.4) 
where the explicit forms of the horizon radii are given by 

r^=(j(l + 2cos(^ + l.)]] , (3.5a) 

r^= (^J(^l + 2cos^-) 1 1 , (3.5b) 
( m ( ^ ^ [6 2 



with r\ + r"^ = —Tq. Here 

6 = arctan 



l + +_,,|| , (3.5c) 



/ 3V3g V4m¥ - 27g^ 
I 2m3Z2 - 27g4 



^In order to define the horizon, the charge q must satisfy a relation q* < 4m^l^ /27. 
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The positions expressed by r_|_ and r_ correspond to the outer and inner horizon, respec- 
tively. It is useful to notice that the charge q can be expressed in terms of 9 and m by 



— ■ sm^ I - 



27 

The outer horizon takes a value in 



—l<rj^< y/ml, 

where the upper bound and the lower bound correspond to the case for q = and the 
extremal case, respectively. 

The temperature is defined from the conical singularity free condition around the 
horizon r_|_, 

where 

a = —Q-, b=—. (3.7) 
zr_|_ 

In the limit q ^ 0, these parameters go to 

l3/2 



a — > 0, 



2mi/4' 

and the temperature becomes 

It might be useful to rewrite the parameters a and b in terms of the temperature and the 
chemical potential, 

b = (^) ^ (3.8) 



1 + Vi + 4(A/r)2 Vvtt; i + ^i + 4(;i/r) 



where we denote fl = fiy^S-KG^ / {3{TTeiy) . 

For later convenience, we shall introduce new dimensionless coordinate u = r'^/r'^ 
which is normalized by the outer horizon. In this coordinate system, the horizon and the 
boundary are located at n = 1 and n = 0, respectively. The background metric ( 3.2a| ) can 
be rewritten as 

with 

f{u) = (1 -ti)(l + n- an^). 

Now let us explain the interpretation of this RN-AdS spacetime in terms of the brane 
setting. We consider Nc D3-branes and treat them as a gravitational background i.e. 
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AdSs X 5^. In order to consider the baryon charge at the boundary theory, we introduce Nj 
D7-branes. The D7-branes are wrapping on S"^ of S^, and we neglect this dependence. 
Here we use the probe approximation for the D7-branes, so that a backreaction from "the 
D7-brane tension" is neglected. Then the effective action becomes that for 5D gauge 
theory. The baryon current corresponds to U{1) gauge field on this flavor D7-branes. The 
5D effective action of the system might be given by the action ( |3.1| ) with 

_ I _ NcNf 

The baryon charge which can be identified by Q and the charge of RN-AdS q are then 
related by ^ 

By using the relation ( |3.11 ), one might understand that AdS leads RN-AdS through a 



backreaction from "the U{1) baryon charge" on D7-branes. However this does not neces- 
sarily imply that Nc and Nf are of the same order of magnitude. The baryon charge is 
provided by open strings between D3- and D7-branes. We can introduce large numbers of 
these strings even for small Nf. Then, the geometry receives the backreaction from the 
charge even though the D7-branes itself is treated as a probe. We here consider the case 
in which Nf is finite while the charge Q is large. This can be understood in the large Nc 
context through the relation ( 3.11| ). The finite contribution of the charge q of RN-AdS 



could be only considerable if the large baryon charge Q{(x y/Nc) is taken in the large Nc 



We will evaluate the quark number susceptibility using the hard wall model [28, 29 1 



and soft wall model pQ] with the RN-AdS background. Note that the Hawking-Page 



type transition in both models is studied in |31]. In the hard wall model to ensure the 
confinement a sharp infrared (IR) cutoff is introduced, while in the soft wall model a dilaton- 
like field is coupled to the 5D bulk action. At high temperature, due to the Hawking-Page 
type transition discussed in |^, AdS black hole background is stable. In this case the 
black hole horizon completely covers up the IR cutoff of the hard wall model, while in the 
soft wall model an IR scale, which is associated with the dilaton-like field, is still visible. So 
the only energy scale in the hard wall model at high temperature is the temperature itself. 
In the soft wall model, however, we have two energy scales, temperature and the IR scale. 
A consequence of this is that the quark number susceptibility in the hard wall model turns 
out to be temperature- independent, while that in the soft wall model shows non-trivial 
temperature dependence at high temperature. The latter is close to the observations made 
in lattice QCD and also in QCD models. 

3.1 Quark number susceptibility in hard wall model 

In this subsection, we shall discuss the quark number susceptibility in the hard wall model 
on the RN-AdS background. The 5D action of the gauge field, which is dual to the 4D 
quark number current j^(x) = q{t,x)'y^q{t,x), is 

1 
47i 



S = -—l d^x^i^^^F-", (3.12) 
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where is the 5D gauge couphng constant. In this work we consider two different values 
of the gauge couphng constant: g'^ = 127r'^ /Nc from the hard wall model |28, 29] and 
g1 = AT^H/{NcNf) from D3/D7. 

To obtain the quark number susceptibility, we use the action ( |3.12| ) with the metric 
and gauge perturbations around the classical background ( p. 2a ) and ( |3.2b| ). Since the 
Green function which provides the quark number susceptibility ( |2.7| ) is the current-current 
correlator Gt tik), we here need to consider the scalar type in the metric perturbation. We 



follow the procedure in |24| and refer to the result 



Gt t(w, k) 



Ik' 



9a 



+ 



2(2 -a)26- 



where 



with 



451(1 + a)62 U2 - 3tj2 ' Dp{uj,k) }' 
Dpioj, k) = 2(2 + a)bk'^ - + a)uj - (2 - afbD^uj'^, 

27(2 - a)aVl + 4a 



(3.13) 



(2-a)4(l + 4a)3/2 



1 



AMl) = -log(2-a)-l„g(l--^^ 



/ /2 + 5a 



4(l + 4a)^/2(l + a)3log(2-a) 
4(2 - 2a + 41a2)(l + a)2Ki(l)|, 
2a 



By using the formula (|2.7D , we obtain 



2c/f 62 y2 + a 
where a and b are given through the definition (|3.8|), 



(3.14) 



b- 



ttT\ 1 + 1 1 + 



Wl 



3(iVc55) 




2 - AnTb. 



(3.15) 



In the charge free case /U = 0, our result agrees with that in [13]. 

Before going on further with the quark number susceptibility, we briefly discuss the 



Hawking-Page type transition. As in ]31], there exists Hawking-Page type transition in 
the hard wall model and soft wall model. At low temperature in confined phase, thermal 
AdS dominates the partition function, while at high temperature in deconfined phase, AdS 
black hole geometry dominates. Therefore, the quark number susceptibility is described 
by the AdS black hole background at high temperature and by the thermal AdS at low 
temperature. The critical temperature for deconfinement is ~ 120 MeV in the hard wall 
model and Tc ~ 190 MeV in the soft wall model ]31]. A similar critical temperature was 
estimated as Tc ~ 210 MeV by using the soft metric model [^]. Note that the value of the 
critical temperature depends on the number of quark flavors, quark mass and quark number 
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density []33| , To obtain the critical temperature in the present case, we have to do the 
Hawking-Page type analysis with charged thermal AdS and RN-AdS backgrounds. Since 
the charged thermal AdS background has not been discovered, we could not precisely fix 
the critical temperature for the deconfinement transition. For simplicity, we assume that 



the critical temperature is around the value obtained in [31]: Tc = 100 MeV in the hard 
wall model and = 200 MeV in the soft wall model. 

The result in ( |3.14| ) is shown in Figure ||, where the gauge coupling of D3/D7 has been 
used. We confirmed that the result with the gauge coupling from the hard wall model 
shows a similar behavior. Below Tc we adopt the quark number susceptibility calculated 



in the thermal AdS background |13|. As well known, the thermal AdS would not exhibit 
the temperature dependence. Therefore, the quark number susceptibility would be a tem- 
perature independent constant, and it turns out zero fl^. For high temperature regime, 
we use different values of the chemical potential, /i = 0.0 ~ 0.8 GeV. Approaching Tc from 
high temperature regime, we observe in Figure || that x<j/^^ shows a blow-up behavior as 
we increase /i, which may indicate the existence of the CEP in the QCD phase diagram. 
Compared to lattice QCD, our study has an advantage that we do not need to do a Taylor 
expansion with respect to the chemical potential. However, our study has a limitation that 
we cannot study the temperature dependence of the quark number susceptibility below the 
critical temperature, which is a generic problem in a model study based on the AdS/CFT 
due to large Nc nature. Also, the RN-AdS is describing a QCD-like or QGP-like system. 
For instance the relation between the chemical potential and its conjugate charge (or num- 
ber) density is given by Q ~ M^^j see (|3.3D . In a realistic system like QGP, the relation in 



general takes the following form Q ~ afi"^ + bfiT'^, where a and b are constants. 
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3.2 Quark number susceptibility in soft wall model 



We shall work in the soft wall model |30|. We consider the following action 

S = --^ / d'xV^ e-*F^„F™-, (3.16) 

with the dilaton-like field ^{x) together with U{1) gauge field Am{x). We make use of the 
RN-AdS geometry (|3.2a|) and ( ^.2b ) as the background. 

In ||3^, the dilaton like field is given by <I>(x) = /^c/r^. Rewriting it in terms of u, we 
have 

<^{x) = cu, (3.17) 
where c = l^c/r"^. We shall work in Au{x) = gauge and use the Fourier decomposition 

A,it, z,u)=j^^ e--*+^^-yl^(A;, u), (3.18) 

where we choose the momenta which are along the z-direction. Equations of motion with 
respect to At{u) and Au{u) are given as, 

= 4' + ^^—^K - {k^At + kLoA,) , (3.19a) 
e * uj \ J 

= uA't + kfA'^, (3.19b) 



where the prime stands for the derivative with respect to u. An equation of motion for A^ (u) 
can be derived from ( p. 19a ) and ( |3.19b ). For Ax{u) and Ay{u), one can obtain decoupled 
second order ordinary differential equations. Since we are interested in the time-time 
component of the retarded Green function to calculate the quark number susceptibility, we 
will not consider Ax{u) and Ay{u) hereafter. 

3.2.1 Solution for At 

From the equations ( t3.19a ) and ( 3.19b| ), we obtain an equation of motion 

= A+ -^Z^At + (^(a; - k f) - (3.20) 



Since the equation ( 3.20 ) is an ordinary second order differential equation for A')-(u) with a 
regular singularity at the horizon u = 1, we first write the solution as A'f.{u) = (1 — uYF[u) 
where F{u) might be a regular function at the horizon. The parameter u could be fixed as 

v = -i—^, (3.21 
47rr' ^ ' 

by imposing the incoming wave condition. 

Now we solve the equation of motion in the hydrodynamic regime i.e. small oj and k 
compared with the temperature T. In order to do the perturbative analysis it might be 
useful to factorize the dilaton factor from F{u), 

F{u) = e*(")F(u), (3.22) 
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so that the equation ( |3.2C1| ) can be expressed as a simpler form 



= ( e*u/F' 

26 
2^a 



+100- e*ii(l + u — av?)F' + iuj-^ (e*n(l + u — au^) ) F 



+ 



uH^ e^ 
(2 — a)^ 1 + ti — an- 



■((2-a)2 + (l-a)(3 — a)u 



+ (1 - 4a + a^)^^ - a(2 - a)^^ + j F 

-kH^e'^F. (3.23) 

The function F{u) is now expanded as 

F{u) = Fo{u) + ujF^{u) + k'^Fk2{u) + 0{uj'^ ,ujk'^), (3.24) 

where the functions Fq(u), Fi^(u) and Fi^2 (u) are determined by imposing suitable boundary 
conditions. The solution can be obtained recursively^. The result is as follows^: 



Fo{u) = C, (const.) 

F^{u) = iCbeH Ei{-cu) + Ki{u) - Ki{0) 



(3.25a) 



2- a 



Ei{c{l-u))-E,{c)-logil-u)) }, 



(3.25b) 



Fk2{u) 



2„c 



E,{-cu) + Ki{u) - Ki{0) 



2- a 

(1 + a)e 



E^icil - u)) - Ei{c) - log(l - u) 



2(2 - a)Vl + 4a 

(1 - a)e 



log 



/ l-2au\ 



1 + 



1 



log 



/ l-2au\ 

1 



-e log u + 



2(2 - a) 



log(l + u — au^) 



(3.25c) 



with 



Ki(u) 



2(2 -a)^/ 



^ , |e-^(^+v™) f (1 + a) - (1 - a)VTT^) 



^The derivation of the solutions is given in Appendix^. 
^The function Ei{x) is an exponential integral 

/■°° -t 
Ei{x) = - dt . 



This can be expanded as follows: 



where 7 is the Euler constant. 



Ei{x) =7 + loga:: + a; + C'(a::^), 
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xEi(—{l + VI + 4a - 2au 



- Vl + 4a - 2au) 

Let us consider the integration constant C. This could be estimated in terms of the 
boundary values of the fields 

lim At{u) = A°t, lim A^{u) = A^,. 

Using the equation of motion ( t3.19a| ), a relation 

lim (ufiA'{ - <I>'A[)) = (^2^0 + ojkA^j (3.26) 

should hold. Therefore we may fix the constant C as 

b(k^A^t+iokA^) 
C = ^ T — . (3.27) 

e'=fia;--(l-e-")A;2' 



c 

One can see the existence of the hydrodynamic pole in the complex w-plane. 

Near the boundary the obtained solution leads a relation between the radial derivative 
of the fields and its boundary values 

A'^{u) = b^i^k^A^t + ujkA^^ logn 

+ — b{k^ A'l + LokA^)i +O{ojk'^,k'^) \+0{u). (3.28) 

(iio - -{I - e-'')^ 

One can also obtain a similar relation for A'.^{u) through the equation ( |3.19b| ). 
3.2.2 Retarded Green functions 

We now proceed to evaluate the Minkowski correlators. An on-shell action can be obtained 
from ( 3.16| ), 

u=l 

. (3.29) 

By using the relation ( |2.4D and the definition (p.5|), we can read off the correlators in the 
hydrodynamic approximation, 

J k^ 

2glbe^ {ioj-Dk^)' 
I ujk 
2^e^ {iLO-Dk'^) 
J 

2glbe^ {ioj-Dk^)' 



So[A] = / 7^ e-*(«) (Ati-k,u)A[ik,u) - f{u)A,{-k,u)A',{k,u) 



Aglb^J (27r) 



Gt tiuJ,k) = TT^^T^ -—2^, (3.30a) 



Gt ziuj,k) = -TT^^T^ -72V' (3-30b) 



G, ,{u,k) = TT^^T^ — (3-30c) 
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where we have introduced the fohowing local counter term in the boundary to remove the 
logarithmic singularity: 



^5 

The constant D is the diffusion constant 



5ct = ^\ogej d'x^f^)F^,F^''. 



D = -{l-e-'). (3.31) 

c 



Using the correlator (3.30a), we can obtain the quark number susceptibility Xq iii terms of 



the temperature and the chemical potential through the definition (p. 7]), 
In the charge free case it reduces to 

which is different from [13|. We confirmed that ( |3.33| ) is correct starting from the 5D 
AdS-Schwarzschild background. In terms of c(= cr'^/l^), the susceptibility ( p. 32]) is given 

by 

2lc 

In the soft wall model [^], c is fixed by the rho meson mass. In the present case we 
cannot use hadronic observables such as masses or couplings to fix the constant since we 
are working in the black hole phase, where light mesons are to be melted away. Here we 



take another route to fix it. We compare our C2 defined below, equation (|3.35D , with that 



from lattice QCD [14| at T = Tc and choose c to reproduce the lattice result. 



X,/r2 = ^2n(2n - l)c2„(/i/r)2(""i). (3.35) 

n 

The quark number susceptibility with the chemical potential is shown in Figure |2|, 
where the gauge coupling from D3/D7 has been used. For high temperature regime, we 
use different values of the chemical potential, fi = 0.0 ~ 1.6 GeV. Again we find that the 
quark number susceptibility shows a blow-up behavior as we lower the temperature to the 
critical temperature, thereby indicating the existence of the CEP in the (^, T) plane QCD 
phase diagram. 

3.3 Implication in QCD phase diagram 

Here we discuss our results in the light of QCD phase diagram. As well known, the nature 
of the QCD transition does depend on the number of quark flavors and the value of the 
quark mass. For pure SU{3) gauge theory with no quarks, it is the first order. In the 
case of two massless and one massive quarks, the transition is the second order at zero or 
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Figure 2: Xg/iNcNfT'^) in the soft wall model for varying /i(GeV) with iVc = 3 and Nj = 2. 



small quark chemical potentials, and it becomes the first order as we increase the chemical 
potential. The point where the second order transition becomes the first order is called 
tricritical point. With physical quark masses of up, down, and strange, the second order 
at zero or low chemical potential becomes the crossover, and the tricritical point turns into 
the critical end point, see ||T^ for reviews. 

Now what can we say about the QCD phase structure based on our study? Our 
approach is based on AdS/CFT, and so it inherently probes QCD (or QCD-like theory) at 
large Nc- The transition suggested by our study is manifestly first order as shown in the 
figures, |l| and|2|, where the quark number susceptibility shows a discontinuous jump at the 
transition point. This discontinuity is there since we use thermal AdS at low temperature 
and RN-AdS at high temperature. This first order nature of the transition would be the 
end of story as long as we don't consider 1/Nc corrections. Though we may not be able to 
assemble all the leading 1/Nc corrections in a consistent way, we discuss some speculation 
on what could come out with those corrections. For zero chemical potential, as shown in 
Figs. H and we have two lines at low and high temperature regimes. There is a hope that 
with 1/Nc corrections those two straight lines are connected with no discontinuity at zero 
or small chemical potentials since the discontinuous gap is not much big. If this happens, 
the transition could be the second order or the crossover. When we increase the chemical 
potential, this possibility does not seem plausible since the quark number susceptibility 
in high temperature regime blows up near the transition point and also the transition 
point will move to the left, i.e. the transition temperature decreases with the chemical 



potential [33, 34 1. This means that the quark number susceptibility at low temperature, 
which is zero with no 1/Nc corrections, should increase very rapidly to realize smooth 
change at the transition temperature. With 1/Nc corrections, this will not be possible. 
Therefore at high chemical potential, the first order nature of the transition will persist 
regardless of the presence of 1/Nc corrections, while for small chemical potentials it may 
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change to the second order or the crossover due to the l/N^ corrections. However, we 
emphasize here again that this is just a speculation. 

In short, the transition from our study is the first order, and it may, however, become 
the second order or the crossover with \/Nc corrections. 



4. Quark number susceptibility under magnetic field 

In this section, we study the quark number susceptibihty with an external magnetic field 
turned on. The basic motivation is due to the observation that the constant magnetic field 
enhances the dynamical chiral symmetry breaking {qq) ~ | eB \ and generates the dynamical 
quark mass rriq^^ = f{\eB\) |25|. Using an effective low energy QCD model, linear sigma 
model with quarks, the authors of [^5| showed that with increasing magnetic field the QCD 
transition changes from the crossover to the first order. This implies that there exists the 
CEP in {B,T) plane as we raise the magnetic field. Therefore, we would expect that the 
external magnetic field affects the behavior of the quark number susceptibility at finite 
temperature. For instance, the presence of the CEP means a peak in the value of quark 
number susceptibility as we increases the magnetic field. So the main motivation of our 
study with the magnetic field is to observe if the peak appears with the magnetic field. In 



addition, a recent study |26| shows that sufficiently large magnetic fields are likely created 
in relativistic heavy ion collisions, and so our study may be tested in experiments. 

To calculate the quark number susceptibility at finite temperature, we consider here 
the non-extremal AdSs x S^, 

ds^ = ^li^ ( _ f{u){dtf + {dxf) + ^!^(dn)2 + fd^l (4.1) 



where /(n) = \ — and T is the temperature. The gauge field comes from the probe 
D7-brane whose action reads 

Sd7 = -NfT^ j d^x e-'t>yJ-dei{gMN + 2-Ka'TMN): (4.2) 

where T^ = l/((27r)'^/f ) and e*^ = Qs- Qmn{x) is the induced metric which we consider as 
the trivial one. Here the external magnetic field B enters the action as p6|] 

J'mn = Fi?]v + Fmn, i^i? = B. (4.3) 
Then, wrapping the D7-brane on S'^ and taking ^ 



Emn = Qmn + F^^\, (4.4) 



5 



we get the 5D action 

= -^/d'^ dnV-det(i?5) {E^^T\E^^f^FmnFik, (4.5) 

''One may be tempted to use the hard wall or soft wall model for simplicity. In this case, however, the 
magnetic field does not affect the equations of motion for gauge fields. 

^Since (-E^^)"'" is not symmetric, there exist additional terms in the action, which are the powers of 
(E^^)"^" Fnm.- However, our choice of 4-momentum, fc^ = {uj,0,0, k), makes those terms vanish. 
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where = y/4:TT'^l/{NcNf) and we have defined an inverse of 5D part of Emn i-e. 
{E^^)'^'-E^in = E^rii{E^^Y^ = 5™. As we did in the previous section, we take the gauge 
-^u{x) = and the same Fourier decomposition as ( |3.18| ). The resulting equations of motion 
with respect to At{u) and Au{u) lead 

= X^,\XtuA!t)' - {k^At + kioA,), (4.6a) 
= uA't + kfA'„ (4.6b) 

where we have defined 



Xt, = V-det{E,){E^Y{E^'y^, X,, = V-det{E,){E^Y{E^T- (4-7) 

Differentiating the equation ( [4.6a| ) with respect to u and using the equation ( 4.6b ), we 
obtain 

= [xi,\XtuA',)')' - {k' - io^r')K (4.8) 
Now using the hydrodynamic expansion 

A'tiu) = {l-uY[FQ{u)+ujF^{u) + k^F^,2{u) + 0{uj'^,ujk^)y 

where v = —iuj/{ATTT) as the incoming wave condition, we get the equations of motion for 
Fq{u) and Ff.2{u), respectively, 

= {Xi,\XtuFo)'y, (4.9a) 
= {Xi,\XtuFk2)'y - Fo. (4.9b) 

In ( |4.9aD , Xf„Fo(ti) should be a constant (= Co) to avoid the singularity at u = 1 due to 
XfJ^{u) — )• as n — )• 1. As a result, we obtain 

Fo = CoXrJ. (4.10) 



Using this solution, the equation ([4.9bD is recasted as 

Xi^\XtuFk2y = CoS{u) with S{u)= j^du'Xt^\u'). (4.11) 

Then, we insert the above solutions into ( [4.6a| ) to obtain 

k'At + kujA, = X;,\XtuK)' 

= Xi,\Xtuk^F,2y + 0{uj) 

= k^CoS{u) + 0{oj). (4.12) 

Since we will take uj = 0, this equation determines At{u). Note that At{u) obtained from 
this procedure is the zeroth order term of the series solution, since we are substituting 
A'l{u) evaluated to 0{uj,k'^). In fact, we should integrate the equation ( 4.11|) once more 
to get the solution of 0{k'^), which cannot be done analytically. However, to get the 
susceptibility, only the zeroth order solution is needed since A; — >• limit of Green function 
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with a; = implies the contribution of the zeroth order only. Thus, from ( 4.121) , the zeroth 
order solution is now 



At{u) = CoS{u) = A^l 



5(0) 



(4.13) 



where ^° = Af,{u)\u=o. 

As a result, we get the retarded Green function and the quark number susceptibility 



Xg = - limRefCt t(w = 0, A;)) = ^[5(0)] 
Then, using the explicit form of Xtu{u) from ([4.7|), we obtain 



-1 



(4.14) 



9l 



9l 



° dn 





-/dn 



47ra'S 



arcsinh 



27ra'5 

(^^IT)2 



-1 



(4.15) 



When = 0, we observe that Xq/T"^ ~ (constant), which agrees with the result of 
the hard wall model at /i = 0. This should be so since the hard wall model action used 
is nothing but the leading term of Dirac-Born-Infeld action for the D7 probe brane after 
perturbative expansions of the action in terms of a' . 

The plot of Xg/^^ for varying B{= 27ra'B) is given in Figure where we take B = 
(0, 1,2,3,4) from bottom to top. We find that the quark number susceptibility increases 
rapidly with increasing B as we lower T from high temperature regime. Note that since 
D3/D7 model does not exhibit the confined phase, we plot our quark number susceptibility 
in entire finite temperature regime. This observation itself is very new, though the blow- 
up behavior at low temperature could be expected by the study of [35|. By studying 



the magnetic field dependence of a modified potential, the authors of |35] showed that 
with increasing magnetic field the QCD transition changes from the crossover to the first 
order, which implies the existence of the CEP, and so the diverging behavior of the quark 
number susceptibility. In this sense, our study support indirectly the result obtained in 
|35]. But, our study made in this section would not address the QCD transition itself due 
to the absence of the confined phase in D3/D7 model. To improve this defect and to see 
if the blow-up behavior is universal regardless of the gravity background, we may consider 
D4/D6 or D4/D8 model, which is relegated to a future study including the effect of a finite 
quark mass. Note that the finite quark mass seems soften the peak in the quark number 
susceptibility [^] . 

Before closing this section, we compare Xq/T"^ iii D3/D7 and D4/D8 |^8| at finite 
temperature with no magnetic field. We prepare a basic setup to calculate the quark 
number susceptibility with an external magnetic field in D4/D8 in Appendix]^. In p9|| , 
from the point of view of the chiral symmetry breaking, an external magnetic field was 
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Figure 3: Xg/iNcNfT'^) for varying B = ina' B with ^ Z and Nf ^ 2. Here B = (0, 1,2,3,4) 
from bottom to top. 

considered in D4/D8. Although we do not provide the result in Appendix B, we can see 
the behavior resulting from D4/D8 system by look through the B = case. In that case, 
Chern-Simons terms do not contribute and we can easily get the result 



-du 



u~ 



(4.16) 



This result is, however, different from what we obtained in D3/D7 system. The different 
T dependence of Xq come from the different T dependence of the horizons and from the 
different exponents of u in the integrands. The dimensionful parameter / compensates the 
different powers of T. On the baryonic density and susceptibilities in D4/D8 model, we 



refer to [|0|. 



5. Summary 

We studied the quark number susceptibility with the finite quark chemical potential or 
under the external magnetic field at high temperature. We first considered the hard wall 
model in the RN-AdS background and observed that as we lower the temperature starting 
from a high temperature, the quark number susceptibility shows a peak with increasing /x. 
When ^ = 0, however, the quark number susceptibility is independent of the temperature 
both at low and high temperatures. To improve this at least in high temperature regime, 
we move to the soft wall model and found that at high temperature the quark number 
susceptibility increases with the temperature for fi = 0, and observed a similar blow-up 
behavior as we lowered the temperature with increasing chemical potential. This peak 
may imply the existence of the CEP in QCD phase diagram on (fi, T) plane. We discussed 



implication of our results to the QCD phase diagram in 3^ , where we concluded that the 
transition from our study is the first order, and it may, however, become the second order 
or the crossover with 1/A^c corrections. 
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We then calculated the quark number susceptibility under the external magnetic field. 
As we raise the magnitude of the magnetic field, we observed a similar rising-up behavior 
as we lowered the temperature with increasing magnetic field. It will be interesting if this 
behavior in {B,T) plane is to be confirmed or disconfirmed by lattice QCD or any other 
studies. 

Finally, we compared our D3/D7 model study with D4/D8 at high temperature to find 

X, ~ D3/D7, 
X, ~ T\ D4/D8. 

We observed that the T dependence of Xq is sensitive to the background geometry. In QGP, 
Xq/T'^ saturates to a constant value, ideal gas limit, at sufficiently high temperature. It is 
interesting to see that D3/D7 model shows this feature. This might imply that the blow-up 
behavior in {fi,T) or {B,T) plane may vary with the gravity background. Therefore it will 
be interesting to see if the blow-up behavior is universal. On top of this, finite quark mass 
effect on the quark number susceptibility is also interesting to investigate. These issues 
will be addressed in a future study. 

Before closing the paper, we remark some limitations of our study. In our study, we 
can calculate the susceptibility with arbitrary values of the chemical potential. However, 
we are not able to study the temperature dependence of the quark number susceptibility in 
confined phase due to the Hawking-Page type transition. Therefore, we evaluate the quark 
number susceptibility only in high temperature regime, deconfined phase. In addition, our 
study based on AdS/CFT is inherently suffering from corrections. Unfortunately, a 
systematic way of collecting all those corrections has not been established. 
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A. Perturbative solution of At in the soft wall model 

Plugging the expansion ( |3.24| ) into the equation ( 3.23| ), one can read off the one for Fq{u), 



= [e^ufF^j . (A.l) 

A general solution is given by 

Fo{u) = Co + Do\Ei{-cu) - ^Ei{c{l - u)) + i^i(n)|, (A.2) 
I- 2 — a J 

where Cq and Dq are integration constants. Imposing the regular condition at the horizon, 
the solution should be 

Fo(n) = Co = C. (const.) (A.3) 
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The constant C could be estimated later. 

By using this solution, one can get an equation for F^(u), 

0= (e*u/F^y + iw^^(e*n(l + u-atx2)yC. (A.4) 

A general solution is 

F^{u) = Ci + Di{Ei{-cu) + ifi(n)} 

-^^{DiEi{c{l - u)) - iCbe'logil - n)}. (A.5) 
Removing the singularity at the horizon, the integration constant Di should be 

In order to fix the another constant Ci, it might be convenient to impose a condition at 
the boundary, 

F^{u) - Ei{-cu) lim ' 



u->o \ Ei{—cu) J 



0, 



u=Q 



so that finite terms at the boundary could be removed. The final form of the solution is 



F^{u) = iCbe^i Ei{-cu) + Ki{u) - i^i(O) 



■{Ei{c{l - u)) - Ei{c) - log(l -u))\. 



2-a 

A differential equation for Fj^i{u) is 

0= (eVi^fe2)'-C6V. 
A general solution can be obtained as 

Ffc2(u) = C2 + D2\^Ei{-cu) + Ki{u)^ 

--—{D2Ei{c{^ - u)) + log(l - u)} 



+ 



l + a 



2(2 - a)c Vl + 4a 



log 



/ 1 - 2au \ 



1 — 2au 



+2(2 - a) log u-{l-a) log(l + u- av?) \ , 



(A.6) 



(A.7) 



(A.8) 



and the constant D2 can be fixed as 



Do = - 



Cb^e' 
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The remaining constant C2 might be determined to ehminate finite contributions at the 
boundary. The solution then becomes 



Ffc2(n) 



C6^e' 



2„c 



E,{-cu) + Kiiu) - Ki{Q) 



2- a 

(l + a)e 



E,{c{l - u)) - Ei{c) - log(l - u] 



log 



/ l-2au\ 



2(2 - a)Vl + 4a 

fl - a)e-^ 



1 + 



1 



log 



-e log li + 



2(2 - a) 



log(l + u — au^ 



/ 1 - 2an \ 

1 



(A.9) 



B. D4/D8 brane system with the external magnetic field 

In this section we shall discuss an effect of a constant magnetic field to the quark number 
susceptibility. We here consider the D4/D8 brane system with the constant magnetic filed. 

We first introduce the bulk background geometry of Nc D4-branes in type IIA super- 
string theory with the compactification on a circle. There exist two different phases i.e. 
confined and deconfined phases. We here take the deconfined phase. The background is 
then given by 

d.^ = {^"\ - fmdtf + (df)^ + (dx4)^) + (§)'^'(7^ + u'^^i) > (B.la) 

= 9s [^-^j , (B.lb) 

where the </)(x) is the dilaton field. dO^ is the metric of the four sphere and R is the 
curvature radius of the background which is expressed in terms of the string coupling Qs 
and the string length Is = ^fa' , 

R^ = ng^NJl 

The function f{U) is given by 

/(C/) = l-(^y, (B.2) 
and the temperature can be read off as 



3 a 



1/2 



T 



4tT i?3/2 • 



(B.3) 



Following Sakai and Sugimoto |38|, we introduce the probe D8-brane which sits in the 



transverse direction to the compactified one X4. In the deconfined phase where the X4-U 
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subspace forms a cylinder, the D8-brane might be a straight hne which simply follows the 
geodesic from U = Ut to infinity, 

dsls = [^y\ - fmdtf + (dxf) + {§y\^ + U'dni^ . (B.4) 

The action for the D8-brane consists of the sum of the DBI and the Chern-Simons 
actions. The DBI action is given by 

Sdbi = -Ts j d\ e-^ V- det(W + 2T:a'FMN), (B.5) 



where Tg = l/((27r)^Zg) is the D8-brane tension and Qmn{x) is the induced metric (|B.4| ). 
We put the constant magnetic field B as the background of the U{1) gauge field in the 
D8-brane and consider small fiuctuations, 

J'mn = Ff^^ + FMN with Fi°) = S. (B.6) 

It might be useful to collect the background fields as 

Emn = Qmn + 2vra'F];]y . (B.7) 

Integrating over the four-sphere, we can then obtain the following action for the fiuctuations 
from the DBI action 

S5D = J d'xdUy^-det{E,) i^-j {E^'r\E^' f"F^nFik, (B.8) 

where the indices m and n run through t, 1, 2, 3, U and an inverse of 5D part of Emn has 
been defined i.e. (E^^)"^^ E^in = E^niiE^^Y"^ = ■ We set the four sphere components of 
the gauge fields to be zero. The 5D Chern-Simons action arises after an integration of the 
RR four form over the four sphere on the D8-brane, 

Sc^ = -i-^\AKFKF. (B.9) 



^487r2^ 

As we did in the main part of the paper, we work on the gauge Au{x) = and use the 
same Fourier decomposition as (3.18). Equations of motion for At{U) and Az{U) are then 
given by 

° ^ "{ ( © ^^^^^0 " © 'tM ^ ""^^'^ } " ^^'^^""^ 

= «|(^(^)5(C/)/(f/K^ +(|) ^j^^[uj'A,+ukA,)^-ipBA[, (B.lOb) 
° " " (^) ^^^^ ("^^^ ^ ^fiU)A'.)^ - i/3B (ojA, + kAt) , (B.lOc) 



where 



9{U) = \l{^y + {27ra'B)\ {B.U) 
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The constants a and /3 are defined by 

" 247r3a" ^ 87r2 ' 
The equations ( [B.lOaD and (|B.10cD imply ( |B.10b| ). 

In order to solve the set of equations, it is standard to introduce the master variable, 

Z{JJ)=u:A^{U) + kAt(lJ). (B.12) 

The the master equation becomes the following form: 



V g{U) f{U){u^-ky{U))) \R) p{U) 

_.PB fuy ' ukf'jU) 
' a {rJ f{U)g{U){oj^-k^f{U)) 

It might be convenient to introduce a dimensionless radial coordinate u = Ut/U. The 
locations of the horizon and the boundary correspond to u = 1 and u = 0, respectively. By 
using this coordinate, the master equation ( B.13D becomes 

^[u^ g{u) ^ /(n)(a;2 - k^f{u)) ) ^ [AttT J up{u) 

+ '3 'B ^^f'i"^) ^ 

uf{u)g{u){uj^ - k^f{u)) 

+{37Ta'Bf—-^-—-, (B.14) 
u^fiu)g^{u) 

with 



w 



where the prime now implies the derivative with respect to u. We can impose the incoming 
wave condition at the horizon, 

Z{u) = {l-u)-'4^F{u), (B.15) 

where the function F{u) should be regular at the horizon. The master equation then 
becomes that for the function F{u), 

\u giu) f{u){uj^-k^f{u)) 2ttT1-uJ 
f w 1 / 1 1 g'(u) oj'^f'iu) 



AttT \ — u\l — u u g{u) f{u){oj'^ — k'^f{u)) 
uj{u)g{u)(uj^ - k^f{u)) u^f{u)g^{u) 
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Multiplying the factor (w^ — k'^f[u)) to the equation above, we could apply the hydrody- 
namics approximation. The function F{u) can be expanded as 

F{u) = FQ{u)+ujF^{u) + kFk(u)+0{iJ,e,ujk). (B.17) 

The order 0(ujk) in the expansion of the master equation, we can fix the function -Fo('") as 

Foiu) = 0, (B.18) 

Equations for Fk{u) and Fi^{u) can be read off from the 0{uj'^k) and 0(w/c^), respectively, 

= (uf{u)g{u)F'f,{u))' + i3Tra'Bf'{u)FUu) + {STra'Bf^^, (B.19) 

= (u9{u)FM)' + i37ra'B(-^]' Fk{u) + i37ra'Bf ^"^"^ (B.20) 
^ ^ \fwJ uj(u)g(u) 
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